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Abstract 

We develop the basis of the two dimensional generalized quan- 
tum statistical systems by using results on r-generalized Fibonacci 
sequences. According to the spin value s of the 2d-quasiparticles, 
we distinguish four classes of quantum statistical systems indexed 
by s = 0,1/2 : mod(l), s = 1/M : mod(l), s = n/M : mod(l) 
and < s < 1 : mod(l). For quantum gases of quasiparticles with 
s = 1/M : mod(l), M > 2,, we show that the statistical weights den- 
sities p M are given by the integer hierarchies of Fibonacci sequences. 
This is a remarkable result which envelopes naturally the Fermi and 
Bose statistics and may be thought of as an alternative way to the 
Haldane interpolating statistical method. 
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1 Introduction 



Quantum elementary excitations carrying generalized values of the spin, to 
which we shall refer here after as quasi-particles, are involved in the descrip- 
tion of many aspects of the critical behaviours of physical systems in two and 
two plus one dimensions [1,2]. In effective field models of fractional quan- 
tum Hall (FQH) liquids for instance [3], the elementary excitations carry 
fractional values of the spin and fractional electric charge Q. In this kind of 
systems, the spin angle 6 is a topological quantity which is generated by a 
Bohm Aharonov effect and eventually by a Berry phase for systems with non 
zero space curvature geometries. For FQH hierarchical states of generic level 
n, typically described by a U(l) n CS gauge theory with a matrix coupling 
K of integer entries K; L j 1 < i,j < n and integers k, the spin angle ^ , 
proportional to s = Yhj=i h{K~ x )ijlj, is not necessary half integer. The spin 
s takes rational values and in general real ones. 

From the mathematical point of view, the parameter s may, roughly 
speaking, be thought of as the highest weight of the 5*0(2) representations. 
As such there is no restriction on it and, a priori, can take any value. This 
parameter, together with the scaling dimension 5 were at the basis of the 
spectacular developments made over the two last decades in 2d critical sys- 
tems [4,5,6]. 

In 2d conformal field theory (CFT), objects carrying generalized values of 
the spin are elegantly described by world sheet (WS) vertex field operators 
type V a ((f)) = exp(ia(p); where describes one or a collection of 2d free 
scalars and a is a parameter related to the square root of the spin s of V a (4>) 
[7]. Quantum states with general values of the spin are also involved in the 
description of the continuous spectrum of the Liouville theory [8], in the 
integrable thermal deformation of minimal models [9] including the 7Ln Potts 
ones [10] and those related with U q (sl 2 ) quantum group representations [11] 
especially the periodic ones [12]. 

In statistical mechanics of many bodies, quasi-particles carrying gener- 
alized spins have been subject to some interest in the past. In [13]; these 
quasi-particles have been interpreted as virtual particles interpolating be- 
tween bosons and fermions. This is explicitly seen on the Haldane statistical 
weight W(a) giving the number of accessible states of a gas of iV quasi- 
particles occupying a group of G states: 

W(a) = (1) 
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In this equation, a is an interpolating parameter lying between zero and one 
describing respectively bosons and fermions. Later on we will give the exact 
relation between a and the spin s. 

Starting from this statistical weight formula, Wu has studied in [14] the 
thermodynamic properties of gases of quasi-particles with exotic spins and 
derived the basic equations of the general form of the distribution and the 
grand partition function. These equations are exactly solved for bosons, 
fermions and for the so called semions. Some of the results of the Wu study 
were rederived in [15,16] by using a completely different manner based on an 
extension of the Pauli exclusion principle [17]. The key idea of this method 
may be summarized as follows. Considering a gas of quasi-particles of spin 
1/M mod 1, M integer greater than one (M > 2) and assuming that the 
quantum statistics of quasi-particles of spin 1/M are governed by a gener- 
alization of the Pauli exclusion principle according to which no more than 
(M — 1) quasi-particles can live altogether on the same quantum state. This 
principle extends the usual Pauli principle for fermions, associated here to 
the special case M — 2. In [15,16], see also [18], it has been shown that 
one can solve exactly several statistical mechanics quantities for this kind of 
quantum gases. The obtained results are in complete agreement with the 
Wu analysis [14]. 

In this paper we reconsider the study of [15,16] by examining the basis 
of such kind of quantum statistical systems using results on combinatorial 
aspects of generalized r-Fibonacci sequences. Instead of the interpolation 
formula of Haldane eq(2), we propose a new statistical weight p(N + G, M), 
given by the following combinatorial expression of M-generalized Fibonacci 
sequences: 

P(n, M) = V ( * + ?t"" + '7- ) W ■ • ■ , (2) 

g +2g l +-+Mg M _ 1 =n-M 9o-9l- ' ' ' 9M-1- 

where the p^s, < Pi < 1, are interpreted as occupation probabilities. To 
have an idea on how this relation applies in practice, it is instructive to 
consider the pure case where the p^s are set to one and so the above formula 
reduces to the well known M-generalized Fibonacci numbers namely: 

p(n,M) = £ G/o + gi + --- + gM-i)! - (3) 

go+2 gi +-+Mg M -i=n-M 90-9l- ' ' ' 9M-1- 

Moreover setting g + g± + ■ • • + gu-i = G or equivalently g\ + 2g 2 + • • • + 
(M — l)gM-i = N, one obtains the generalized weight density of a gas of 
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quasi-particles of spin 1/M, whose leading case M — 2, gives just the usual 
fermions statistical weight W(l) = p(N + G, 2)\ go+gi=G = N uQi N y ■ 

This paper is organized as follows. In section 2, we give the set up of 
quantum statistics of quasi-particles carrying general values of the spin and 
fix terminologies. In section 3, we review some general features of the Haldane 
Statistics and compare it to our approach. In section 4, we make an explicit 
study of the quantum gas of quasi-particles of spin s = 1/3 either by using 
standard techniques or in terms of M — 3-generalized Fibonacci numbers. In 
section 5 we review the M-generalized Fibonacci sequences and exhibit their 
connection with fractional statistics. In section 6 we make a discussion and 
give our conclusion. 



2 Generalized QSM 

Statistics is a quantum mechanical feature of identical particles which governs 
the quantum behaviour of a collection of iV particles in the large limit of N. 
Usually one distinguishes two kinds of quantum statistical systems: bosons 
carrying integer values of the spin (s = modi) and fermions carrying half 
odd integers spins (s — 1/2 modi). The thermodynamic statistical properties 
of bosons and fermions are described by two different quantum statistical 
distributions namely the Bose and Fermi densities. In this section we want 
to discuss the general set up of the quantum statistics of systems of quasi- 
particles carrying general values of the spin. Then we give some results 
regarding the quantum statistics of quasi-particles of spin 1/M, M > 2. 

To start recall that in two dimensions, there is no restriction on the values 
the spin s; it can take any value; i.e: < s < 1 mod 1. Thus, according to the 
values of the spin s, the quasi-particles describing the quantum excitations 
in 2d systems may be classified in different kinds and so it is natural to 
conjecture the existence of various types of quantum statistical systems in 
two dimensions. These are given by the following classification: 

• QSM; quantum statistical systems with quasi-particles carrying spins 
s = 0,1/2 (modi). 

• FQSM; quantum statistical systems with quasi-particles carrying frac- 
tional values of s; i.e:s = 1/M : mod(l), where M > 2 
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• RQSM; quantum statistical systems with quasi-particles carrying spins 
s type: s = n/M : mod(l), where M > 2 and M ^ k.n with k an 
integer. 

• GQSM; quantum statistical systems with quasi-particles carrying spins 
s, with < s < 1 : mod(l). 

This classification, according to the values of the spin s, obeys naturally 
to the following embedding: 

QSM c FQSM c RQSM c GQSM. (4) 

It is not difficult to see that QSM, describing bosons ( s = mod (1)) and 
fermions ( s = 1/2 mod (1)), is recovered as the boundary limit of FQSM by 
taking respectively M = +oo and M — 2. Similarly, FQSM may be viewed 
as a particular case of RQSM by setting n — 1. Finally, putting s = n/M in 
GQSM by restricting to quasi-particles with rational values of the spin, one 
obtains RQSM. 

From the above classification, we expect that several features of QSM 
may be extended iteratively to the other statistical systems. For example 
the well known fermionic number operator N = [(— ) F ] extends naturally as 

< ip | N | ip >=< if; | expi(2nF) \ ip > (5) 
= expi(2ns), (6) 

where | ip > is the wave function of a given quasi-particle state of spin s. 
The eigenvalues of the operator F are effectively given by the spin of the 
| ip > eigenfuntions; namely s = : mod(l) for bosons, s = 1/2 : mod(l) for 
fermions, s = 1/M : mod(l) for FQSM, s = n/M : mod(l) for RQSM and 
< s < 1 : mod(l) for GQSM. 

An other example of features of QSM which extends naturally to the other 
quantum systems of eq(4) is the Pauli principle. This exclusion principle 
may be generalized to all other quantum statistical systems. This is easily 
seen by computing the trace of fermionic number operator expi(2nF) on 
representation spaces involving all possible values of the spin; in particular 
for the supermultiplets built of states (s = 0, s = 1/2) for QSM and their 
generalizations (s = j/M, : < j < M - 1) for FQSM. 

Tr : expi(2nF) = 1-1=0, for QSM, (7) 
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M ~ X 27T7 

Tr : expi{2ixF) = ]T [expi( — j)} = 0, for FQSM, (8) 

i=o 

Such relation extends to RQSM and GSSM as well. For the supersym- 
metric multiplets (s = jn/M, : < j < M — 1) of the RQSM systems; we 
have 

M ~ l 2-Kin 

Tr : expi(2nF) = ]T [expi( — = 0, for RQSM, (9) 

3=0 M 

while for the multiplets (s — f3, : < f3 < 1) of GQSM systems, we have the 
following formula: 

Tr : expi(2TrF) = f 1 ds[expi(27rs)} = 0, for GQSM. (10) 

Having classified the various kinds of quantum statistical systems, we turn 
now to describe briefly methods used in studying gases of quasi-particles. 
There are several ways of doing. The first way, due to Haldane, is based on 
the combinatorial formula eq(l) which turns out to be the simplest and the 
natural weight one may conjecture. This approach uses linear interpolating 
relations between the Bosons and the Fermions and concerns quasiparticles 
of spin s lying between 1/2 and 1. Haldane quasiparticles obey a generalized 
exclusion Pauli principle although a clear statement of this principle for the 
continuous values of s is still lacking. A second approach which was suggested 
in [19], follows the same line of reasoning as of Haldane; but instead of 
looking for interpolating statistical weight, one considers interpolating wave 
functions. In other words if we denote by (f> and ip respectively the Bose and 
Fermi wave functions of a system of N particles, then the wave function \ °f 
gases of quasiparticles of spin s is, roughly speaking, given by: 

X = (j)cos(2ns) + i/jsin(27rs); < s < 1/2. (11) 

However, a detailed analysis of this eq shows that for the case of quasi- 
particles of spin s = 1/M one discovers exactly the quasi-determinant of 
[15,16] describing the wave function of the FQSM gas of spin s = 1/M (mod 
1); with M > 2. 

A third way, which was developed concerns gases of quasiparticles of spin 
s = 1/M (mod 1), is based on the idea of imposing from the beginning 
the generalized Pauli exclusion principle according to which no more than 
(M — 1) quasiparticles of spin 1/M; M > 2 can live altogether on the same 
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quantum state. In this approach, one obtains explicit relations for FQSM 
extending the Bose and Fermi analogue; in particular the following quantum 
distributions 

^ ^ * M ^ M Mo x 

< n k >= y—^ - j3^m, (12) 

where we have set £ = e~^ €k ~^a,nd M > 2. This eq offers an alternative to 
the standard quantum q-uon one given by < n(q) >= Yr^;where — 1 < q < 1 
and where for q — 1, — 1, we get respectively the Bose, Fermi and Boltzman 
distributions. . 

In the next section, we review the main lines of the Haldane statistics 
eq(l).Then, we develop further the alternative we considered in [15,16]. Us- 
ing group theoretical methods and combinatorial analysis we show that the 
rigorous results we obtained are not fortuitous but are due to the existence 
of a link between FQSM systems and a well known subject in mathematical 
physics namely Fibonacci sequences hierarchies. To that purpose, we shall 
proceed as follows: 

(a) first review the basis of the Haldane approach; this method is based on 
a linear interpolating formula between bosons and fermions. The simplicity 
of this way is due to the fact that the weight density is due to the remarkable 
properties of the special function j3(x,y). Non linear interpolations between 
bosons and fermions were also considered in [19]; they involve hypergeometric 
functions J- {a, b; a, f3) and. 

(b) Develop a new setting using rigorous analysis involving group theo- 
retical techniques and combinatorial formulas and show that the obtained 
results are indeed linked to generalized Fibonacci sequences. 



3 Haldane Statistics 

We start this section by describing some general features of the Haldane 
statistics. Although exact solutions are still missing, this statistics exhibits 
many remarkable features which make it very special. Haldane statistics 
based on eq(l), corresponds in fact to a linear interpolating statistics between 
bosons and fermions. This linearity is a useful mathematical feature which 
make the interpolating formula natural and simple to handle. Indeed, setting 

x = x(a) = N , 
y = y(a) = G-aN + (a-l), 1 j 



8 



in eq(l), one sees that it reduces to the inverse of the well known di-gamma 
special function (3(x,y); namely: 

W(x,y) = fc±fl, (14) 
x\y\ 

W(x,y).p(x,y) = 1. (15) 

Note that y(a) is a linear interpolation between the Bose point ys and the 
Fermi one yp respectively given by: 

y B = y(l) = G-N ( , 

y F = y(0) = G-l. 

The interpolation eq(13) as well as the statistics of the boundary points 
eq(16), shows that the Haldane interpolating parameter a and the spin s are 
related in the large N limit as 

Ml -2s) , . 1 
a = N _ 2 ~ (! - 2s), <s<-. (17) 

For s = 0, and s — |, the parameter ct is equal to one or zero respectively. 
For the special case a = |, corresponding to semions of spin s — |, the 
Haldane statistics can be solved exactly; one finds the following quantum 
distribution for semions is given by: 

<„,(„= 1/2) >= ;?T -^— __ (18) 

For general values of a, the situation is however complicated because one has 
to solve non linear eqs of the type: 

1 

< n i > 



£ = ™(0 Q [1 + ™(0] 1_Q (19) 
where £ = e /3 ^ e_ ' 1 - ) , e is the energy and \x is the chemical potential. 



4 Gas of Quasiparticles of Spin 1/3 

The gas of identical quasi-particles of spin 1/3 is the simplest system coming 
after the quantum Fermi gas. This is a quantum system of FQSM, which 
may be studied by using two routes: 
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(i) using the Haldane interpolating statistical weight W(a) eq(l) and 
fixing the parameter as a = 3 (^_ 2 ) ' i- e s = 1/3. 

(ii) solving the constraint equations imposed by FQSM according to which 
no more than (M — 1) quasi-particles of spin s — 1/M can live altogether on 
the same quantum state. 



4.1 Interpolating Statistics 

Following Haldane and using eq(17), the statistical weight W s= i/ 3 of a gas of 
iV quasi-particles of s = 1/3 can be shown to be given by: 

r[ 3(iV-2) G !-(* -D(2iV-6) ]r(iv) 

VV s=i / 3 - 3(^-2)0-2^-2^+6 1 • y* U > 

1 i 3(7V-2) \ 

This eq may be rewritten into an equivalent form as: 

W 9=1/3 = W(N,y), (21) 

where now y = 3( - jV 1 _ 2 ) — 2)G — N 2 + 6 — 2N\. Having these relations 
at hand one can derive the thermodynamic statistical mechanical features of 
this gas of quasi-particles. For example, the mean occupation number < > 
is obtained by solving eqs(19). For the case of a gas with pure occupation 
probabilities, the partition function z(£) reads, at the statistical equilibrium, 
as: 

*(0 = i+e+e 2 - (22) 

Note that eq(22) lies between with the partition functionz Fermi(0 — 1 + £, 
and z Bo se{0 = 1 + Y^=\i n of individual fermions and bosons of quantum 
Fermi and Bose gases respectively. ZFermi(0 is given by just the sum of 
two terms as a consequence of the Pauli exclusion principle contrary to the 
bosonic gas where no exclusion rule exist. Therefore, the polynomial eq(22) 
of order 2 reflects clearly that quasiparticles of spin s = 1/3 obeys a gener- 
alization exclusion statistics principle according to which no more than two 
particles of spin s = 1/3 can live together on the same quantum state of the 
gas. This selection rule is the simplest exclusion statistics coming after the 
Pauli exclusion formula obeyed by fermions and constitutes the leading term 
of a more general hierarchical formula namely: 

^) = l + £ + ... + £ M - 1 . (23) 
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This relation extends eq(22) for the case of quasi-particles of spins -p, M > 2. 
In what follows we will show that such relation may be also derived by using 
an other alternative using combinatorial analysis. 



4.2 Generalized Exclusion Statistics 

We start by examining the special exclusion statistics obeyed by quasi- 
particles with spin s = 1/3 by using group theoretical methods and combina- 
torial analysis. Then we work out the properties of the formula of the statis- 
tical weight at equilibrium one gets, in particular its relation to 3-generalized 
Fibonacci numbers. To do so, consider a system of N quasi-particles of spin 
s — 1/M with M — p + 1, occupying G quantum states. On each state of 
the G ones, there can live altogether at most p quasi-particles as required 
by the generalized exclusion statistics for particles with spin s = l/(p + 1). 
To get the total number of accessible state of the p = 2 system, we proceed 
as follows: First divide the set of the allowed G quantum states into three 
kinds of sub-sets go, g\ and g% consisting respectively of Bo, B\ and B2 boxes 
which can contain zero, one and two quasi-particles of spin s = 1/3 of the 
system with occupation probabilities po, p\ and P2 with < pi < 1. Thus, 
we have the following equalities: 

£0 + £i+#2 = G, (24) 
9l + 2g 2 = N. (25) 

Since the B , B\ and B 2 boxes are indistinguishable, it follows that for a 
given configuration of the quantum system; the number of accessible states 
is no longer given by Glp^p^p 92 but rather by: 

G ' ! iPfvfpf- (26) 



For configurations associated with pure occupation probabilities; i.e po = 
Pi = P2 = 1, the above relation reduces to the well known — — r factor. 
Taking into account all possible partitions of the number N , the total number 
D 3 (N, G) of accessible states of the quantum gas occupying G quantum states 
reads then as: 

D 3 (N,G)= £ , ' ^>fpf, (27) 

9U92>0,gi+2g 2 =N 9lW{^ ~ 9l ~ 9lV- 
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where we have used eq(25). Furthermore solving gi in terms of iV and g 2 by 
help of eq(26), one may rewrite this equation into following equivalent form. 



W/2] r , 

= g gt(N _ 2gm + g _ N / M'P?, (28) 

where we have set g = g2 and where [N/2] is the integer part of N/2. Note 
in passing that this eq has some remarkable features; one of which is that 
D 3 (G,N) contains the fermion statistical density D 2 (G,N) as the leading 
term of the sum over the possible partitions of the B 2 boxes. In other words; 
D 3 (N,G) may be written, for pure probabilities, as 

MN > G) = N\(G-N)\ + g g\(G-2g)\(G + g-N)\> (29) 

where the first term of the right hand side (rhs) of is just the statistical weight 
D 2 (N, G) of N fermions occupying G states. The remaining ([iV/2] — 1) extra 
terms of eq(30) describe then a deviation above the D 2 (G, N) density and so 
it may be interpreted as the effects of quasi-particles of spin s = 1/3. Put 
differently, we have 

D 3 (G, N) = D 2 (G, N) + A 32 (G, N), (30) 

where A 32 = D 3 (G, N) — D 2 (G, N), given by the second term of the rhs of 
eq(32), describes the manifestation of FQSM. 



5 Fibonacci Sequences Setting 

A careful observation of the previous combinatorial analysis shows that it has 
much to do with the combinatorial realization of the generalized Fibonacci 
sequences [20], see also [2 1,22]. To exhibit this connection explicitly, let us 
first give some elements on the combinatorial aspects of these sequences, in 
particular for the Fibonacci sequences of order M=2 and 3. To start let do, 0,1 
be some fixed real numbers with a\ 7^ and consider the sequence {V^}^, 
defined by V n — p{n + 1, 2) with, 

p(n, 2)= £ ^{^' agV, , :: for every :: n > 3, (31) 

S0+2ffi=n-2 90-9l- 
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and p(2, 2) = 1, p(n, 2) = for n = 0, : 1. Using this eq it is not difficult to 
verify that V = 0, Vi = 1 and 

K+i = a V n + aiK-i, for every : n > 1, (32) 

implying in turns that {V^}£^ is nothing else a weighted Fibonacci sequence 
with initial conditions specified by the numbers Vq,Vi and weights ao,a\. 
Note that for ao = a± — 1, we recover the classical Fibonacci numbers, 

p{ n , 2 ) = 51 , :: for ev ery " n > 3. (33) 

go+2gi=n-2 #0-<7l- 

Observe also that up on imposing the extra constraint g + g\ = G, the V n is 
intimately related to the fermions number density T> 2 (G, N). Indeed putting 
the condition go + g\ = G back into eq(36), one discovers that sum over g\ 
and go with g + 2g± = n — 2 is completely fixed and son = iV + G' + 2as 
well as: 

V 2 (G, N)=p(N + G + 2, 2) | fl0+fll=G , (34) 

showing that the fermion density V 2 (G, N) is given by a restricted Fibonacci 
numbers sequence of rank two. An inspection of this result shows that it is 
not a fortuitous feature and may be used to describe not only fermions but 
also statistical weights of quantum gases of fractional spin s = 1/M. This 
observation is supported by several arguments, above all the fact that Fi- 
bonacci sequences admit integer hierarchies in perfect coherence with FQSM 
systems. To illustrate the idea, we consider in what follows the M = 3 case; 
i.e the case of a quantum gas of quasi-particles of spin 1/3. We begin how- 
ever by giving some generalities on 3-generalized Fibonacci sequences, then 
we derive the relation linking D 3 and p(n, 3). Later on we give the general 
result. Let a , ai, a 2 be some fixed real numbers with a 2 ^ 0. If a 2 = but 
ai 7^ 0, one falls in the previous case. A 3-generalized Fibonacci sequences 
{Kijn^o is roughly speaking defined by: 

K+i = a V n + aiK-i + a 2 K-2, for every : n > 2, (35) 

where Vo, Vi, V 2 and a , a±, a 2 are respectively the initial conditions and the 
weights of {V^}£^. Eq(38) extends the usual Fibonacci sequences {V^}^, 
obeying eq(44), and which may also baptized as the rank 2 Fibonacci se- 
quences; the leading element of the Fibonacci hierarchy. . Recall by the 



13 



way that r-generalized Fibonacci sequences, which have been studied by us- 
ing various method and techniques, go beyond eqs(44) and (47), and play a 
fundamental role in many fields of applied sciences and engineering [ ]. In 
practice there are two ways to handle these sequences: the combinatorial re- 
alization which we consider hereafter and the Binet representation. For the 
third order of the hierarchy we are discussing, the combinatorial realization 
is generated by: 

p(n, 3) = Yl — T^ «"f, , " for every :: n > 3, 

ffo+2si+3ff 2 =n-3 90-91-92- 

(36) 

where we have set p(3, 3) = 1 and p(j, 3) = for < j < 2. To check 
that the above eq is indeed a 3-generalized Fibonacci sequence, it is enough 
to take V n = p{n + 1,3). Straightforward calculations show that we have 
Vq = V\ = and V 2 = 1 as well as: 

V n+ i = a V n + a{V n -i + a 2 K-2, for every : n > 2. 

Note in passing that for a = a x = a 2 = 1, the sequence {V n }n=o gives the well 
known 3-generalized Fibonacci numbers, very familiar in the mathematical 
literature [M],[K]. In this case, eq (45) becomes: 

P(n, 3) = E ^±gL±fH , :: for every :: „ > 3, (37) 

go+2<?i+3<? 2 =n-3 90-91-92- 

with p(3, 3) = 1 and p(j, 3) = for < j < 2. The beauty of this relation 
is that it is linked to the density of the quantum gas of quasi-particles of 
spin s = 1/3 we studied earlier. Indeed, if instead of the two constraint eqs 
(17-18), we consider rather the following relaxed condition 

g + 2 9l + 3g 2 = G + N, (38) 

obtained by adding eq(17) and eq(18), then the density D 3 (G,N), given by 
eqs(28-29), is translated exactly to eq(46). Put differently D 3 (G,N) can 
be derived from the 3-generalized Fibonacci numbers by imposing either 
9i + 2(fe — N or g + gi + g 2 = G. In the formulation based on the Fibonacci 
sequences, the number of accessible states eq(19) is now replaced by: 

( 90 + 9l + 92)\ p9 9 9 
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Taking into account of all possible configurations of the number N + G parti- 
tions in eq(41), the new total number p 3 (N, G) = p(G + N — 3, 3) of accessible 
states is: 

Ps(N,G)= £ (g °7\ + f 2) Po>rpf, (40) 

g Q +2 gi +3g 2 =G+N 00-01 -02- 

which sould be compared with eq(39). Introducing the Dirac delta function 
S(x), we can write, 

D 3 (G, N) = 8(g + gi + g 2 - G)p(G + N — 3,3). (41) 

Substituting the delta function and solving gi in terms of N and g%, we can 
rewrite eq(20) into the following equivalent form: 

[N/2] 

where we have set g = gi and where [N/2] is the integer part of N/2. Having 
this relation at hand, we can write down all thermodynamics properties of 
the gas of quasi-particles at equilibrium. 



6 More on FQSM 

In this section we want to extend some of the previous results to quantum 
systems of iV quasi-particles of spin 1/M; M > 2 using the combinatorial 
M-generalized Fibonacci sequences. We start by generalizing the useful prop- 
erties on the combinatorial and the Binet expressions for the M-generalized 
Fibonacci sequences. 

Let do, di, om-i be (M — 1) arbitrary real numbers with Om-i 7^ and 
M > 2. An M-generalized Fibonacci sequence {V n }n=o is defined as follows: 

V n+ i = a Q V n + aiK-i + ••• + a M -iV n - M +i, for : n > M — 1 (43) 

where Vq, Vi, Vm-i and ao, ai, bm-i are the initial conditions the se- 
quence. In case where = 1 for all i's and V — Vi — ... — Vm- 2 = 
and Vm-i = 1, we recover the M-generalized Fibonacci numbers while for 
a M-i = and om_2 ^ 0, one gets .M — 1-generalized Fibonacci sequences 
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to which one can apply the same analysis. The combinatorial expression of 
{Ki}n^o> f° r n > M, is realized as: 

Pin, M)= £ + f |--- + ^- l)! ag°af • • • , 

g +2g 1 +-+Mg M _ 1 =n-M 90-9l- ' ' ' QM-1- 

(44) 

where we set p(M, M) = 1 and p(n, M) = 0for0<n<M — 1. Starting 
from the following identity 

(A; + --- + A;m-i)! _ ^ (to + • • • k M -i - 1)! 



k \ki\ ■ ■ ■ fc M -i! ~^ fc ! • • • kj-ilikj - l)\k j+ i ■ ■ ■ k M -i ] - 

we derive easily that, 

p(n + 1, M) — a p(n, M) + a lP {n — 1, M) + h a M -ip(n - M + 1, M), 

Thus the sequence {W n }n5) defined by W n = p(n + 1,M), is a sequence 
whose initial conditions are such that W — Wi — ... — W M -2 = and 
Wm-i = 1- More generally, the combinatorial expression of the sequence 
(46) can be given, for every n > M, in terms of p(n, M) (47) as follows. 

V n = A p(n, M) + A lP (n — 1, M) + h A M -ip(n — M + 1), (45) 

where 

A m = a M -iV m + a M _ 2 V m+1 H h a m 1/ M -i; for :: m = 0, 1, • • • , M - 1, 

and where p(M, M) = 1 and p(j, M) = for j < M. 

In the Binet representation of M-generalized Fibonacci sequence, the V^s 
are expressed as polynomials in the eigenvalues of: 

P{X) = X M - a Q X M - 1 a Af _ 2 X - a M -i, (46) 

The point is that to every sequence {V n }^=o of rank M, we associate the 
polynomial P{X) above 

whose roots A , Ai, • • • , \ s have respectively the multiplicities m , mi, ■ • ■ , m s 
The Binet realisation of the general term V n of the sequence reads in terms 
of the A-s as: 

s m s — 1 
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where the f3' l -s are obtained from initial conditions by solving the following 
system of M linear equations 

s m s — l 

E E Pij ni tf = V m ■ for : n = 0, 1, M - 1. (48) 

Such expressions are known in the mathematical literature as the Binet for- 
mula of {V n } n > . For the leading n=2 case where the Xj solutions of the 
characteristic polynomial P(x) = x 2 — a x — a± have are simple to handle, 
we have 

V n = j^BjX], for A^A 2 , (49) 

V n = (Si + B 2 n)A?, for A x = A 2 , (50) 
where the B/s are related to «o and cti as 

2 

E S A H = a «5 n = 0,l(for A^A 2 ) (51) 
j'=i 

(S! + S 2 n)A^ = a n , n = 0,l(for A 1 = A 2 ). (52) 

Such eqs may written in a more explicit way by specifying the explicit 
form of the X'jS in terms of the initial parameters. 

Having given some tools on Fibonacci sequences hierarchies, let us turn 
now to expose briefly the link with FQSM systems. To that purpose, we 
reconsider the gas described in section 4 but now the quasi-particles have 
spin s = l/M with M = P + 1, P > 1. In this case eqs(27-28) extend as 

9o + 9i +92 + •••• + g P = G (53) 
gi + 2g 2 + 3g 3 + .... + Pg P = N (54) 

The density of accessible states of the gas is, 

V P+1 (N,G)= ]T G\[(G-{2G P )\]l(g l \)}- 1 (55) 

9i,92,....,gp>0 k=1 1=1 

Ef =1 j.gj = N 

Eq(37) exhibits several features, some of them extends properties of the 
V 3 (N,G) statistical weight discussed earlier. This concerns for instance the 
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interpretation of V P+ i(N, G) as a deviation above the V P (N,G) weight of 
quasi-particles of spin l/P. Indeed rewriting eq(54) by singling out the sum 
over g P ; i.e, 



[N/P] 



G\ 



V P+1 (N,G)= J2 — E , , ur r<' w 

9p =o9p<- 9i\92\....,9p-i\{G-Y, j= i93)\ 

9l: 92, 9p-l dL U 

Ej=i j-9j = N - Pgp 

(56) 

one sees that the leading term in the sum over g P ; ie for g P — 0, is just the 
weight density T> P {N ) G) of a gas of N quasi-particles of spin l/P (mod P). 
Put differently, eq(55) may usually decomposed as 

V P+1 (N, G) = V P (N, G) + A P+liP ; (57) 

where Ap +1]P measures the derivation from V P+ i(N, G) to V P (N, G) densi- 
ties and reads as 

[N/P] 1 G] 

a p+ i,p= J2 —\ S rm — : — ttt^ — _ 



9p ^9p\ ^....^^q 9p UC v ; l9j )\ 



Xf- 1 L j.g j = N-Pg P 

(58) 

This eq is just a generalization of eqs(32-33); it is valid for any value of 
P greater than three and obeys quite remarkable properties inherited from 
the generalized Fibonacci combinatorial method we have been using. Eq(57) 
may be also thought as a manifestation of tne exclusion statistics in FQSM 
stating that no more than P quasi-particles of spin -p^- can live altogether 
on the same quantum state. It satisfies equally the following splitting 

V P+1 (N, G) = [V 2 (N, G) + X>,(AT, G)} + £ A, +lj ; (59) 

3=3 3=2 

in agreement with the embedding (4). Here the Aj+ij's, given by similar 
relations to (57), are deviations from Dj + i to Dj densities. 

7 Conclusion 

In this paper we have developed a new way to approach generalized quan- 
tum statistical mechanics of gases of quasi-particles. Our method, which 
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is based on group theoretical analysis and combinatorial representation of 
M-generalized Fibonacci sequences, can be viewed as an alternative to the 
Haldane approach. To do so, we have first reviewed the setting of gener- 
alized quantum statistics which we have classified into four classes accord- 
ing to the values of the spin s of the quasi-particles: (1) the usual QSM 
for s = 0,l/2mod(l), (2) FQSM for s = 1/M mod(l), (3) RQSM for 
s = n/M mod{\) and (4) GQSM for < s < 1 : mod{\). This classifi- 
cation obeys the inclusions: QSM C FQSM C RQSM C GQSM. 

Then we have reviewed the interpolating analysis of Haldane based on 
eq(l) and showed that it is equivalent to the method we developed in [14,15] 
using a generalization of the Pauli exclusion rule as a general principle. 

Next motivated by a set of observations; in particular the embedding 
eq(4), we have developed an other way to deal with FQSM gases. Our method 
is based on the following weight density: 

p(n, M) — £ (90 + 91 + - - + ^-0^ . . . ^ 

g +2g 1 +-+Mg M _ 1 =n-M 90-9l- 4 ' ' 9M-1- 

(60) 

where the p^s are dilution probabilities. Up on imposing eqs(24-25) or equiv- 
alently eq(38), we rederive the known results. We believe that our way of 
doing offer a powerful manner to deal with such quantum systems. Our 
motivations for this claim are: 

(1) Eq(60) is obtained by using a rigorous analysis based on a group 
theoretical method and a combinatorial analysis. 

(2) Eq(60) is in agreement with the embedding QSM C FQSM. It re- 
covers standard results of QSM and describes naturally any kind of quantum 
gas of quasi particles of spin s = 1/M, M > 2.. 

(3) Eq(60) is linked to a well known subject in mathematical physics; 
namely the combinatorial aspect of the M-generalized Fibonacci numbers. In 
this setting QSM is described by the standard Fibonacci number sequences 
while FQSM is described by its integer hierarchies. QSM appears then as 
the leading term of the generalized Fibonacci hierarchy. In this context we 
expect that RQSM and GQSM would be also described by means of similar 
analysis. 

(4) Finally eq(60) incorporates in a natural way the dilution property 
of quantum gases carried by the pi parameters. Therefore diluted quantum 
gases of FQSM are described by integer hierarchies of Fibonacci sequences. 
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